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A hybrid ferromagnet/semiconductor device is used to determine a lower bound on the spin 
lifetime for conduction electrons in silicon. We use spin precession to self-consistently measure the 
drift velocity vs. drift field of spin-polarized electrons and use this electronic control to change the 
transit time between electron injection and detection. A measurement of normalized magnetocurrent 
as a function of drift velocity is used with a simple exponential-decay model to argue that the lifetime 
obtained (« 2ns) is artificially lowered by electronic effects and is likely orders of magnitude higher. 



C/5 



C 

O 

o 



> 

oo 
o\ 

CO 

o 
o 



X 



Manipulation of electron spin injection, transport and 
detection in semiconductors promises to form the ba- 
sis of the future's lower-power, higher-speed information- 
processing paradigm. Much work over the past decade 
has elucidated useful techniques in this nascent field 
of spin electronics (or "spintronics" ) but until recently 
they have been applicable only to direct-bandgap ma- 
terials. This is either because they are optically-based 
(electroluminescence, [l], d, H, HJ Faraday effect, d, @ or 
Kerr effect Q), or because they require epitaxial (metal- 
lic or semiconducting) ferromagnet/semiconductor sys- 
tems restricted to the few available. [1, Q This limita- 
tion was particularly unfortunate because it excluded 
the indirect semiconductor Si, the most dominant semi- 
conductor in the modern microelectronics market, and 
a particularly attractive material for spintronics due to 
the expected long spin lifetime. [Tol. ITT1. fl2| For spintron- 
ics to deliver on its promise and have a serious hope 
of integration with conventional information-processing 
devices, considerable progress must be made in study- 
ing and applying spin transport to this material since it 
forms the basis of the CMOS architecture of computing 
circuits. 

Our recent work with spin-dependent hot electron 
transport in ferromagnetic (FM) thin films for spin injec- 
tion and detection has solved this long-standing problem 
and has opened the field of spintronics to include Si, and 
(in principle) many other indirect-bandgap semiconduc- 
tors previously excluded from spin transport study. [14| A 
side-view illustration and a schematic conduction band 
diagram of our device in equilibrium is shown in Figure 
1(a). Hot electrons are emitted from a tunnel-junction 
(biased with voltage Ve) and through a FM thin film to 
vertically inject spin-polarized electrons into the conduc- 
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FIG. 1: (a)Schematic illustration of the Spin transport de- 
vice, including side-view and corresponding conduction band 
diagram, (b) Second collector current, in in-plane magnetic 
field, showing a spin valve effect at Ve = —1.8V, and (c) pre- 
cession and dephasing (Hanle effect) in a perpendicular field 
at Ve = —1.8V and Vci = 1.0V. Note the pronounced shift 
in the first extrema from zero which is caused by in-plane 
stray field; see text. 



tion band of the undoped single-crystal Si layer. This Si 
is the (lOO)-oriented device layer of a float-zone-grown 
silicon-on-insulator (SOI) wafer, so it has very low impu- 
rity and defect concentration. After transport through 
the 10/im-thick Si, the final spin polarization of these 
electrons is detected on the other side of the Si by eject- 
ing them from the conduction band into a second FM thin 
film and measuring the magnitude of the ballistic compo- 
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FIG. 2: (a) Comparison of predicted prependicular magnetic 
field value of first minima (corresponding to an average n pre- 
cession angle) using Drift-Diffusion model (Eq. [TJ and simple 
Drift-only model (Eq. [2j. (b) Spin precession measurements 
in perpendicular magnetic fields at various externally applied 
accelerating voltages Vci- All curves have been shifted to 
40pA at zero field, a value typical for injection conditions 
used (Ve = —1.8V at 85K). The arrow shows the progression 
of the first minima (used to determine drift velocity) to higher 
magnetic field values with higher accelerating voltage Vci.(c) 
Injected current Ici as a function of the externally applied 
accelerating voltage, where the sharp threshold indicates flat- 
band conditions where electric field £ — 0. (d) Combining the 
results gained from (b) and (c) allows us to compare with 77K 
time-of-flight data of charge transport from Ref. [ijj. Tak- 
ing into account in-plane stray field (Eq. [3} gives excellent 
agreement. 



nent of this hot-electron current in a n-Si collector on the 
other side. Typical spin- valve behavior of the second col- 
lector current Ici in an in-plane magnetic field is shown 
in Fig. 1(b), and coherent spin precession and dephas- 
ing in a perpendicular magnetic field (Hanle effect) [15| is 
shown in Fig. 1(c). These measurements unambiguously 
demonstrate coherent spin transport in Si.[l6j 

In Ref. , we used this type of device to estimate a 
lower bound on the conduction electron spin lifetime of 
«1 ns. Here we establish that electronic effects artificially 
suppress this lower bound, allowing the actual intrinsic 
spin lifetime in silicon to be much larger. 

Since we have rectifying Schottky barriers on both 
sides of the Si drift region, applying a voltage across 
it does not induce spurious currents. However, it does 
change the drift electric field and control the drift ve- 
locity and hence transit time from injection to detec- 
tion. The ratio of second collector current difference 
(AIci), normalized by the injected current (Ici) can 
therefore be used as a simple voltage-controlled met- 



ric which is proportional to spin polarization. [14j (The 
normalization is necessary because a slight increase in 
Ici with accelerating voltage increase is caused by en- 
hanced hot-electron collection efficiency at the injection 
interface). [18[ Clearly, if spintronic effects are to blame 
for an increase in Alcijlci with higher drift velocity, 
the functional form (not just percent change) should fit 
to the expected behavior based on exponential spin decay. 
If it does not, then we must conclude that other effects 
arc blame for artificially suppressing the lower bound de- 
rived only from percent change. 

As we originally showed in Ref , the effects of pre- 
cession caused by a perpendicular magnetic field on Ici 
can be used to determine the transit time r from spin 
injector to spin detector. Because Ici is proportional to 
the projection of spin direction on the detector thin film 
magnetization direction, the first minimum in Ici under 
parallel injector and detector magnetization configura- 
tion defines conditions for an average lot = n precession, 
where u> = gfj-sB/h is spin precession angular frequency. 
To make this interpretation more quantitative, we use the 
following expression to model the effects of precession in 
drift (causing drift velocity v) and diffusion (causing spin 
dephasing):[l![20| 

AJ C2 . ~ / P(x = L,t)cos{ut)e- t/T *tdt, (1) 
Jo 

where P(x = L,t) is the solution of the drift-diffusion 
equation describing the distribution function of electrons 
arriving at the analyzer (x = L = 10/im in our device) 
at time t, after spin-polarized injection at t = 0. The 
exponentially decaying component of the integrand takes 
into account a finite spin- lifetime, r s f. 

The magnetic field value of the first minima of this 
expression as a function of drift velocity is shown in Fig. 
[5] (a). Compared with it is the expected magnetic field 
dependence in the absence of diffusion, when r = L/v so 
that 

v = 2L ^ bB * ( 2 ) 
h 

where B n is the perpendicular magnetic field value of 
the first precession minima in a parallel magnetization 
configuration. It is clear that two regimes exist: one at 
higher drift velocity (where drift dominates and the two 
models predict similar precession minima) , and the other 
at low drift velocity (where diffusion is dominant and the 
minima value is close to ~ 50 Oe). The relative percent 
error, given by the dotted line on the right axis, shows 
that there is little difference between the predictions of 
the two models in the drift-dominated regime. Since mea- 
surements we will show here have minima at perpendicu- 
lar fields greater than 100 Oe, we conclude that they are 
firmly in the drift-dominated region. There is therefore 
negligible error in simply finding the position of the min- 
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ima of Hanle spectra and using Eq. [2] to determine the 
drift velocity. 

Typical spin precession measurements at low perpen- 
dicular magnetic fields, revealing the first precession min- 
ima dependence on accelerating voltage Vci, are shown 
in Figure [2] (b). When the applied voltage increases, the 
drift velocity increases with electric field and the min- 
ima move to higher magnetic field. From the position 
of these minima, we calculate the drift velocity from Eq. 
[2] for applied drift-region accelerating voltages between 
-0.3V and +3.0V. 

To correlate these externally applied accelerating volt- 
age to internal electric drift field, we perform a spec- 
troscopy of injected current Ici as a function of acceler- 
ating voltage Vci at constant emitter voltage, as shown 
in Figure (c). The flat Ve=0 V line indicates that in 
the absence of hot electron injection, no current flows due 
to the fully rectifying metal-semiconductor interfaces on 
both sides of the fO/im Si drift region. However, when 
hot electrons are injected over the injector Schottky bar- 
rier with Ve=-1.8 V, a sharp threshold at Vci ~ —0.4V is 
seen. This is where the applied voltage cancels the elec- 
tric field from ohmic voltage drop in the resistive injec- 
tor base thin film, and causes flat-band conditions corre- 
sponding to zero internal electric field. Since the Si drift 
region is undoped, the conduction band is linear so that 
the electric field is constant and linearly related to the 
voltage difference from the flatband point. Therefore, in- 
ternal electric drift field is calculated from the externally 
applied voltage: £ = (OAV + V C i)/10~ 3 cm. 

This internal electric field measurement enables us to 
plot the drift velocity vs. electric drift field. As can be 
seen with the open symbols in Fig. [2] (d), this data taken 
at 85K compares well with time-of-flight measurements 
of charge transport at 77K.[f7| However, at low fields 
there is a significant discrepancy. By noting the shift in 
positions of 7r extrema for parallel and antiparallel mag- 
netization conditions in Fig. Q] (c) (and Fig. 4 (a-c) of 
Ref. [3), we conclude that this discrepancy is caused 
by an in-plane stray field caused by the FM layers. (The 
shift is clearly not due to measurement delay because it 
is less apparent for higher-order precession extrema at 
higher externally- applied perpendicular magnetic fields.) 
When the injector and detector magnetizations are in a 
parallel configuration, stray magnetic flux passes through 
the drift region in-plane, adding in quadrature to the ex- 
ternally applied perpendicular magnetic field, and shift- 
ing the minima to artificially low external magnetic field 
values: 



Binternal — \j ^external + ^ 'stray (3) 

When the magnetizations are antiparallel, the stray 
magnetic field flux is contained by the FM layers and the 
internal magnetic field equals the external magnetic field. 
As can be deduced from Fig. 1 (c), this stray field is on 
the order of B stra y ~ 300Oe, and when it is taken into 
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FIG. 3: (a) Magnetocurrent A/c2 at different electric drift 
fields from in-plane magnetic field measurements. Note ex- 
pected threshold at V/cm. (b) Normalized magnetocurrent 
Alcijlci showing linear behavior with drift velocity, (c) 
Simple exponential-decay model predicting nonlinear behav- 
ior of AIc2/ Ici- The discrepancy between (b) and (c) in- 
dicates that the lifetime lower-bound determined from (b) is 
artificially lowered by electronic effects. 



account using Eq. [3l the correlation between time-of- 
flight and our corrected data is improved as shown with 
closed symbols in Fig. [2](d). 

With these electric field and velocity calibrations, we 
can now attempt to determine AIc2(v) / Ici(v) by mea- 
suring Iqi in external magnetic fields of -45 Oe and +270 
Oe, corresponding to antiparallel and parallel magnetiza- 
tion orientation, respectively (see Fig. 1 (b)). The differ- 
ence between these measurements is shown in Fig. [3] (a) 
for a dense set of electric fields up to «3400 V/cm. This 
spectroscopy displays the same threshold at flatband con- 
ditions that suppresses the injected current Ici (which 
drives Ici) in Fig. [2] (c). The normalized ratio of these 
two data sets is shown in Fig. [3] (b), where the increase 
at higher drift field could be assumed to be the result of 
the correspondingly shorter transit time and hence larger 
final spin polarization according to P = P$exp(—T/T s f). 

To within the approximation set by our signal to noise 
level, A/c2 /Ici appears to be a linear function of drift 
velocity, with the signal magnitude increasing by approx- 
imately 50% from 1 x 10 6 cm/s to 7 x 10 6 cm/s. 

To compare with our simple model, we plot 
exp(—T /r s f), where r = L/v for various values of r s /, 
on the same drift velocity axis in Figure [3] (c). We ig- 
nore the spin injection and detection efficiencies which 
are unaffected by the electric field within the bulk of the 
Si, so the vertical scale is normalized for comparison with 
Fig. [3] (b) . To achieve the observed relative 50% change 
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across the drift velocity range measured in the experi- 
mental data, a spin lifetime of approximately 2 ns must 
be assumed, consistent with our lower bound given in 
Ref. [l4|. However, the nonlinear shape of this modeled 
curve is in striking contrast to the clearly linear behavior 
seen in the experimental data. 

What could be the cause? Because the electric field in 
the drift region changes as the drift velocity (and there- 
fore transit time) change, the spin-orbit interaction which 
transforms an electric field into a magnetic field could 
cause deviations from the simple model. However, spin- 
orbit is very small in Si, and in any case it should act to 
make the nonlinearity in Fig. [3] (b) stronger, not weaker, 
in a real measurement. 

Electronic effects such as emitter Schottky field- 
dependent hot-electron transfer ratio could easily cause 
the observed linear increase in lS.lcij lex-, and we believe 
that these are the most likely source. For instance, as 
shown in Figs. 2 (b) and 2 (c) of Ref [lj], although I C \ 
drives Ic2, the relationship is super- linear, making the 
simple normalization of Iq2 by Ici insufficient to cancel 
electronic effects. 



In conclusion, we have used coherent spin precession 
in a perpendicular magnetic field to determine the drift 
velocity of spin-polarized electrons as a function of ap- 
plied voltage bias in the drift layer of our Si spin trans- 
port device, and used injected current spectroscopy as a 
function of externally applied voltage to determine the 
internal electric field. This allows us to then correlate 
drift velocity and electric field, which agrees well with 
time-of-flight measurements of electrons in similarly un- 
doped Si. [lTj By measuring the normalized magnetocur- 
rent Ici/Ici as a function of drift velocity, and compar- 
ing to a simple exponential spin-decay model, we deter- 
mine a lower bound for spin lifetime of free electrons in Si 
which is limited by electronic effects and is likely orders 
of magnitude higher. [l2j Measurements in devices with 
longer transport lengths will have quadratically higher 
spin-polarized electron transit times. [201] so to probe the 
1-microsec regime at fixed accelerating voltage, the drift 
region should be in the 300-micron range (i.e. the full 
thickness of a Si wafer). 

This work is supported by the Office of Naval Research 
and DARPA/MTO. 
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